Effect of baryon density on parton production, chemical 
equilibration and thermal photon emission from quark gluon 

plasma 



D. Dutta, A. K. Mohanty, K. Kumar and R. K. Choudhury 
Nuclear Physics Division, Bhabha Atomic Research Centre, Mumbai, India- 400085. 

q\~. (February 7, 2008) 

On 
On 



O 

Q 

in 



Abstract 

The effect of baryon density on parton production processes of gg ggg 
and gg ^ qq is studied using full phase space distribution function and also 
with inclusion of quantum statistics i.e. Pauli blocking and Bose enhancement 
■ factors, in the case of both saturated and unsaturated quark gluon plasma. 

m ; 

The rate for the process gg qq is found to be much less as compared to 

Q\ . the most commonly used factorized result obtained on the basis of classical 

ON 1 

approximation. This discrepancy, which is found both at zero as well as at 
Qh| finite baryon densities, however, is not due to the lack of quantum statistics 

Ph. in the classical approximation, rather due to the use of Fermi-Dirac and Bose- 

' Einstein distribution functions for partons instead of Boltzmann distribution 

which is appropriate under such approximation. Interestingly, the rates of 
parton production are found to be insensitive to the baryo-chemical potential 
particularly when the plasma is unsaturated although the process of chemical 
equilibration strongly depends on it. The thermal photon yields, have been 
calculated specifically from unsaturated plasma at finite baryon density. The 
exact results obtained numerically are found to be in close agreement with the 
analytic expression derived using factorized distribution functions appropriate 
for unsaturated plasma. Further, it is shown that in the case of unsaturated 
plasma, the thermal photon production is enhanced with increasing baryon 
density both at fixed temperature and fixed energy density of the quark gluon 
plasma. 
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I. INTRODUCTION 



The future collider experiments with heavy ions will provide an opportunity to investigate 
the physics of strongly interacting matter at extreme energy densities where the formation of 
a new state of matter , the quark gluon plasma (QGP), is expected. With these colliders, it 
may be possible to achieve energy density well above the deconfinement threshold to be able 
to probe the QGP in its asymptotically free 'ideal gas' form. At high energies, the central 
rapidity region may have nearly vanishing baryon density, similar to the early universe, and 
will probably consist of a very dense system of semi-hard partons (mini-jets) which would 
lead to rapid thermalization and extremely high initial temperature M. In the standard 
scenario, the quark gluon plasma formed during the collision is expected to thermalise in 
a typical time scale of (~ 1 fm/c) and the subsequent evolution follows the ideal Bjorken's 
scaling solution |J. Several theoretical calculations based on perturbative QCD approaches 
now suggest that the partonic fluid which is mostly gluonic will attain kinetic equilibrium 
after a proper time To ~ 0.1 — 0.3 fm/c [3]. This is justified given the fact that the rate of 
elastic collisions between quarks and gluons which establish thermal and mechanical equilib- 
rium (referred as kinetic equilibrium corresponding to isotropic momentum distribution) is 
much larger than for inelastic collisions which establish chemical equilibrium. The remaining 
question is whether the high energy density matter remains in QGP phase sufficiently long 
enough to actually reach chemical equilibrium so that it can be identified with the QGP 
phase seen in lattice QCD calculations. 

In the framework of pQCD based models, it is now understood that the densities of 
quarks and anti-quarks may stay well below the gluon density and may not reach up to the 
full equilibration value as required by an ideal chemical mixture of quarks and gluons ||. 
Because of the consumption of energy by additional parton production, the temperature of 
the parton plasma falls down faster than the ideal solution where T 3 r remains constant. 
Essentially, same picture emerges in a somewhat different approach first proposed by Biro 
et al [|] and subsequently by several other authors With a given initial value for 

the energy and number densities, this type of calculations employ ideal fluid dynamics to 
study the subsequent evolution of the kinetically equilibrated QGP phase coupled to rate 
equations to determine the chemical composition. 

As mentioned before, the evolution of the partonic matter toward a fully (chemically) 
equilibrated QGP will be dictated by the parton proliferation mainly through induced ra- 
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diation and gluon fusion Unfortunately, there have been no rigorous thermal field 

theoretical calculations for the rate of quark production since the basic cross section di- 
verges for the mass less quarks. The closest approach to a full calculation has been done 
by using the thermal quark mass and Debye screening mass to regularize the divergence 
HJ9|]. However, as used by Biro et al || and subsequently by many others 0-0, the quark 
production rate gg ^ qq has been factorized using classical approximation where Pauli 
blocking and Bose enhancement factors (effect of quantum statistics) have been eliminated. 
The calculations have also been simplified further by assuming an ideal baryon free plasma. 
There is now enough evidence from theoretical studies that even at RHIC energies, the mid 
rapidity region may not be completely baryon free [JlOj . 



In an earlier work, we had studied chemical equilibration of the QGP at finite baryon 
density where we had considered the processes gg ^ ggg and gg ^ qq as the two 
dominant channels contributing to parton production. It is natural to think that the presence 
of baryon density (//) in quark and anti-quark distribution functions may directly affect the 
process gg # qq, whereas the process gg # ggg may depend on fi only through the Debye 
screening mass used to regularize the divergence. Using modified Fermi-Dirac distribution 
function for quarks (anti-quarks) and a factorized distribution function for gluons valid at 
small baryon density, it was found in the above work that the two rates are not sensitive 
to the baryo-chemical potential particularly when the plasma is unsaturated. Therefore, we 
had used the same factorization (for gg qq) as used by Biro et al in the case of baryon 
free plasma. Although, the parton production rates for unsaturated plasma are not very 
sensitive to baryo-chemical potential, the energy, quark and anti-quark densities (hence the 
rate of chemical equilibration) strongly depend on \x. It was also shown that due to chemical 
equilibration, the baryon rich plasma cools at a slower rate as compared to the baryon free 
plasma. 

In the present work, we have studied the role of baryon density on parton production 
particularly for the process gg # qq in detail using full phase space distribution functions for 
partons and also including Pauli blocking and Bose enhancement factors explicitly. Specif- 
ically, we investigate the effect of quantum statistics particularly in the presence of finite 
baryon density. The present calculations show that the rate for the process gg qq is 
lowered by 20% to 30% as compared to the factorized result of Biro et al as given in Ref. [f|. 
This discrepancy, however, is not due to the lack of quantum statistics, but rather due to 
the use of Fermi-Dirac (FD) and Bose-Einstein (BE) type distribution functions for quarks 
and gluons in Ref. |Q instead of Boltzmann distribution which would have been appropriate 
under classical approximation. With the revised rates, we have also studied the process of 
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chemical equilibration solving a set of rate equations. The quark and anti-quark fugacities 
slow down whereas the gluon fugacity goes up slightly as compared to the previous values. 
However, there is no net effect on the evolution of temperature which remains practically 
unchanged. Further, it is shown that the rates of parton production are insensitive to the 
presence of baryon density particularly when the plasma is unsaturated. The findings of 
the present study are also consistent with our earlier work [[7]. We had earlier studied the 
thermal dilepton yields as a probe of QGP. We investigate here the thermal hard photon 
production as a complementary probe from a chemically non-equilibrated baryon rich quark 
gluon plasma. It is seen that, like the dilepton yields, the photon yield is also enhanced in 
presence of baryon density. 

The paper is organized as follows. In section II, we calculate the parton production 
and chemical equilibration at finite baryon density. In section III, we calculate the thermal 
photon yield from chemically unsaturated plasma both numerically and analytically. The 
results have been summarized in section IV. 



II. PARTON PRODUCTION AND CHEMICAL EQUILIBRATION 

In this section, we calculate the thermal production rates for (i) gluon fusion gg ^ qq 
and (ii) gluon radiation gg # ggg. We also study the process of chemical equilibration by 
solving the ideal fluid hydrodynamics coupled to a set of rate equations. We use the Juttner 
distribution functions for quarks (anti-quarks) and gluons as given by 

qU) + e(w*)/T x q(q)e ±x + eP /T ' h eP /T _ Xg W 

where x = fi/T. The fugacity factor Aj (i = q,q and g) gives the measure of the deviation of 
the distribution functions from the equilibrium values and fi (= /is/3) is the quark-chemical 
potential. The chemical equilibrium is said to be achieved when Aj —>■ 1. As mentioned 
in Ref. 0, the quark (anti-quark) distribution functions could also be written using the 
commonly used definition of fugacities Aq and Aq given by 

X Q(Q) = e±X \{q) ( 2 ) 

so that, there is no need to use the quark chemical potential fi explicitly. However, the 
definition (1) is quite convenient, since at equilibrium where X q =Xq=l, the chemical potential 
associated with Aj vanishes but the baryo-chemical potential still exists. Further, we have 
taken \ q = A 5 at all values of r so that when /i — ► 0, Aq = Aq = X q resulting in a 
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baryon symmetric matter. While we calculate the parton production rates using the above 
distribution functions, we also consider the approximations for quarks (anti-quarks) as 



Jm A 9(? - )e ±- + eP/ T ~ 1 + eP/ T m J ^ 1 ' 

and for gluon distribution function as 

u = vr ( 4 ) 

where /^L = (1 + e p / T ) _1 and /| 9 = (e p / T — l) -1 . The approximation given by Eq.(3) is 
referred as modified Fermi-Dirac (MFD) distribution in Ref. |7j which becomes Fermi-Dirac 
distribution when Xqiq) — > 1. Hence forth, we will refer both Eqs.(3) and (4) as modified 
factorized distributions (MFD) which are different from Boltzmann (BM) type of factorized 
form given by, 

/,(*) = \U)e-^ )/T ; f g = X g e-?l T (5) 



A. Gluon fusion 

We begin by calculating the rate for the process gg qq ||: 

d 1 f ^ f d3 P 2 f d3 Pz f d3 P* to \Wt . \ v 

R ^ =2 J frymJ Wf2Ej j2^2Ej (2^ (2?r) 6 (Pi+ft-^-^) X 

E \Mgg^\ 2 f 9 ( Pl )f 9 (p 2 )[l ~ / 9 ( P3 )][1 - / 9 -( P4 )] (6) 

and 

D 1 /" ^1 /" d3 P? f d3 P* f ^P* to \Wt . \ v 

= 2 7 (2^7 (2^7 (2^7 (2^ (27F) 5 (Pi+ft-^-^) X 

E l^ M | 2 / g (P3)/ $ (p 4 )[i + / s (pi)][i + / S M (7) 

In Eq.(6), the squared matrix element, summed over spin and color, X) \Ai\ 2 is weighted by 
two gluon distribution functions f g for the initial states. The factor [1 — f q ] [1 — / 5 ] indicates 
Pauli blocking for the final states. In the reverse process Eq.(7)), the rate is weighted by 
the distribution functions of quarks and anti-quarks for the initial states and the gluon final 
states each gain an enhancement factor [l + / s ] due to Bose-Einstein statistics. The factor of 
1/2 accounts for the identity of the two gluons. Using the identity that results from Eq.(l), 
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[i - fm\ = ^ 



Q(Q) A 9 



(8) 



and the unitary relation I.M12I 2 = | -A^f 21 1 2 , Eq.(6) and Eq.(7) can be combined to give 



Rg^q Rq^g 



<-gluon 



(9) 



where 



Lgluon 



f d 3 pi 


r d 3 p 2 


f d 3 p 3 1 


1 {2tc) 3 2E 1 J 


(2n) 3 2E 2 J 


(2tt) 3 2E 3 J 



d 3 p 4 



(2vr) 5 (pi +P2-P3- Pi) x 



E \M gg ^\ 2 f 9 ( Pl )f g (p 2 ) f q ( P3 ) fM >< ^ {El+E2) ( 10 ) 

and = T. The above multi-dimensional integral can be reduced further by replacing 

|| = Jd 4 p5(p 2 -m 2 )6(p ) (11) 
and by change of variables to: 



q = Pi +P2, 
P = g(Pi - P2)> 

g' = P3 +P4, 
p' = 2^3-^4)- 



(12) 



One can easily eliminate the energy-momentum-conserving 5 function by carrying out the 
integrals over q' . Using spherical coordinates and defining q along the z-axis where q^ = 
(qo,0,0,q), p^ = (po,psm9,0,pcos9) and p' = (p' ,p' sin (f> sin x,p' sin (f) cos x,p' cos (p), we 
get, 

' — / / d(l i 

-qo/2 



Lgluon 



IdW Jo d% Jo 



rqo/2 rqo/2 roc roc rl rl 

dq I dp I dp' J dp J dp J d(cos6) j d(cos< 



r 2n 

x / dx5 







P 



Pl + 
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P 



go/2 



i2 2 , 

p - m + 
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x 



gg^qqi Jg 



Qo . 

J+ P 0_ 



Qo 

2~ P ° 



fq 



Qo . , 
j+Po_ 



cos 9 
Qo 



Qo Po 



QP 



2" "ft 



8 
0qo 



cosm — 



Qo Po 
QP' 



(13) 



The summed squared matrix element for the process gg — > qq, is given by || (where quarks 
are assumed to be massless) 



ut 2 ,11 t „ 

I (- + - 

3s 2 274 u' 



(14) 



Here 7^ = 2 (AT 2 — 1) and 7^ = 7^ = 2N c Nf are the number of internal degrees of freedom 
for gluons, quarks and anti-quarks respectively. N c = 3 is the number of colors and Nf 
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is the number of massless flavors. Throughout our analysis we have used Nf = 2.0. The 
Mandelstam variables s, t, u used in the above equations are given by 



2 2 

s = % ~ Q , 



t = m 2 



ffl-Ml +2 



u 



2 

(s + t). 



s 1 



4p| 

g 2 



S 1 



smx, 



(15) 

(16) 
(17) 



The matrix element given by Eq. (14) diverges in the limit u, t — > 0. Hence, the medium 
induced effective quark mass plays a crucial role and the divergence of the cross-section can 
be avoided by replacing 'm' with the thermal quark mass defined appropriately for a non 
equilibrium plasma ||TT|Jl2 l 



m = m„ 



roc 

/ dpp[2f g + f q + f q 
Jo 



(18) 



q 3vr 2 

The integrals over p,p',9 and <fi of Eq. (13) have been performed carefully because of 
kinematical constraints to have the limits 

4m 2 



g > 2m, q 2 < (q 2 - Am 2 ), pi < Q 



2 < 2 <f 

4' Po -T 



1 - 



(19) 



Thus we evaluate all but five of the integrals in Eq.(13) trivially || obtaining 
1 



Lgluon 



16(27r)" J2m 



(gg-4m 2 )2 rp* 

dq I dq \ dp 



dp' 



x I "dxT, I^WVjf +Po] /,[f -Po] /j| -Pol e" qo (20) 

where the limits of integration are 



P* 
P* 



4m 2 , i 

v l )2. 

2 y s ' 



Following Ref. 0], we write Eq.(9) in a convenient form 

AnA/ 



Rg^q — Rq^g ~ n g R 2 {l 



\2 



(21) 
(22) 

(23) 



where n„ is the gluon density and R 2 is given by 



R-2 



-L gluon 

AqAq^ 



(24) 
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Fig. 1 

FIG. 1. (a)The quark production rate R2/T as a function of X g (X q = Xq = X g ) with Juttner, 
Boltzmann and MFD distribution function, for x = 0. The solid dot curve is the approximate 
formula used by Biro et al.[4].The dashed dot curve is calculated rate with Boltzmann distribution 
function in the incoming channel and factorised distribution otherwise. (b)Same as (a) but with 



An 



Xg/5). (c) The quark production rate R2/T as a function of x for different values of 



X g where (X q = Xq = X g /5). 



We compute R 2 numerically using Juttner, MFD and BM distribution functions and also 
compare the result with the approximation obtained by Biro et al 
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1 

0.24/V, a 2 s X g T ln{—) (25) 
a s X 



where a s (= 0.3) is the strong coupling constant. Note that the factor A = (X g + cosha;A ? /2) 
which arises due to the thermal quark mass @ is slightly different from what is used in Ref. 
I 

Figs. 1(a) and 1(b) show R2/T as a function of X g obtained numerically using Eq.(24) 
with Juttner, MFD and BM distribution functions at T =0.57 GeV. Initially, we consider 
a baryon free plasma (x = 0). In Fig. 1(a), we assume the gluon and quark (anti-quark) 
content of the plasma is equal (X q = X g ), whereas in Fig. 1(b), the plasma is assumed 
to be more gluon rich (say X q = X g /5). As in Fig. 1(a), the R2/T values obtained using 
Juttner, MFD and BM distribution functions do not have significant difference at all values 
of A 5 . At small X g , R2/T obtained using MFD differs slightly which ultimately merges with 
that obtained using Juttner distribution function, as both Juttner and MFD distributions 
become equal when X q = X g = 1.0 at x = 0. However, when the plasma is gluon rich as in 
Fig. 1(b), R2/T values obtained with above three distributions do not differ much say up 
to X g ~ 0.5, beyond which the MFD and BM results show significant deviation. Recall that 
for X q = Xg/5, MFD and Juttner distributions (for q and q) are not equal when X g — > 1 even 
at x — 0. We can conclude from these results that for unsaturated plasma, the R2 values 
obtained using Juttner distribution do not deviate significantly from the results that can 
be obtained with the factorized distributions like MFD or BM functions. However, these 
results are quite less (20% to 30 %) than what is found on the basis of the so called classical 
approximation i.e. using Eq.(25) given by Biro et al which are also shown in Figs. 1(a) 
and 1(b) for comparison. It is natural to think that the discrepancy found in Biro et al's 
estimation could be due to the classical aspect of the approximation which we investigate 
below. 

In the classical approximation, Eq.(10) can be written as 

Igluon = ~J ^5 J J^[<Tg 9 -+qqVl2\fg{Pl)fg{P2) (26) 

which represents the cross section for the process gg — > qq folded with the distributions for 
the initial particles. The cross section a gg ^ qq - is given by 

The classical approximation assumes Boltzmann distribution function for quark, anti- 
quark and gluon and also eliminates the Pauli blocking and Bose enhancement factors in 
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the final states. The Eq.(26) can be factorized as I g i UO n = n g a 2 /2 where n g is the equilibrium 
gluon density and a 2 is a velocity weighted cross section. Defining R 2 = n g a 2 /2, Biro et al 
estimate a 2 , n g and also thermal quark mass using factorized distribution functions ( same 
as Eqs. (3) and (4) with x=0 which become Fermi-Dirac and Bose-Einstein distributions in 
the limit A, — ► 1.0). Within the present formalism, we can simulate the results of Biro et al 
by replacing f q and f q in Eq.(20) with BM distribution (Eq. 5) and f g with Bose-Einstein 
type (Eq. 4). Also we estimate thermal quark mass using MFD type distributions as in Ref. 
ID (i.e. using Eqs. (3) and (4) at x—0). The simulated result is also shown in Figs. 1(a) and 
1(b) for comparison (see solid curve with dots). The results are in good agreement with that 
of Biro et al's estimation. Therefore, the deviation found using Eq.(25) is due to the use of 
FD and BE type of distribution functions instead of using BM functions which would have 
been appropriate under classical approximation. In fact, our results obtained from Eq.(20) 
using BM approximation is in a way classical. Further it may be mentioned here that the 
first term in the RHS of Eq.(9) arises due to the factors like Pauli blocking (1 — f q ) and Bose 
enhancement (1 + f g ). Only the integral I g i uon is evaluated using various approximations. 
Therefore, even though we use the classical approximation like BM distribution, the final 
results still include the quantum effect. Similarly, Biro et al have estimated I g i uon classically, 
but their final factorization does include quantum effect (see appendix B of Ref. |7| for 
detail). In other words, Biro et al approach will give correct estimation if BM approximation 
is used to estimate R 2 . The above discussions confine only to baryon free plasma. Similar 
discrepancy is found at finite baryon density also. 

Now we investigate the effect of baryon density on R 2 . Fig. 1(c) shows R 2 /T as a 
function of x obtained numerically using actual (Juttner) distribution function at different 
values of fugacities where X q = X g /5. It is interesting to note that x dependence of R 2 is 
rather weak particularly when the plasma is unsaturated. The nearly x independence is 
also evident from Eq. (1). For small values of X q , the contribution from the factor X q e 
in the denominator is not very significant if x is small. The x dependence of the quark 
and anti-quark distribution functions mainly arises due to the e ±x factor in the numerator. 
Since these exponential factors get cancelled in the product, f q f q (hence R 2 ) will have weak 
x dependence at small baryon density if the plasma is highly unsaturated. 

B. Gluon multiplication 

The rate R 3 for the process gg ^ ggg depends on the differential radiative cross section 

m 
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da 3 



da9J dn ° -e(x f - c ^p-)e(VS-k ± coshy) 



(28) 



dq\ dy d 2 k± dq\ dy d 2 k± v 1 k± 

where the first step function gives approximate LPM suppression factor and the second step 
function accounts for energy conservation, s = 18 T 2 is the average squared center-of-mass 
energy of two gluons in the thermal gas. Here k± denotes the transverse momentum, y is the 
longitudinal rapidity of the radiated gluon and q± denotes the momentum transfer in the 
elastic collision. The regularized gluon density distribution induced by a single scattering is 



dn ( , 



C A ot s 



dy d 2 k± ix 2 fci[(k_L — q±) 2 + m 2 D ] ' 
Similarly the regularized small angle gg scattering cross section is 

da 9 J 9 27ia 2 



dq]_ 4 (q 2 ± + m 2 D ) 2 
where Debye screening mass mo given by 

3g 2 roc 



ni 



D 



ir 2 Jo 



dp p [2 f g + N f (f q + /,)] 



where g 2 = 4na s . The mean free path Xf for elastic scattering is then |5],|7J 



-7ia s n g 



rn 



D 



Integrating part analytically 

i-2-k I 



'1 _L_ 2 m D\ 
v 1 t* q t 2 ' 



2tt 







(kj_ - q±) 2 + m 2 D /( fc 2 +q 2 ±+ m 2 D y _ Aq 2 ±k 2 



and defining R 3 = ^n g cr 3 , we can evaluate 



R 3 /T=^X 2 f n g I{X g ,X q ,x) 



where I(X g , X q , x) is a function of X g , X q and x 



l(Xg,X q ,x) 



dx 



s/4m 2 D 



dz- 



o 



1 + * 



COsh (y/x) 



X 



+ 



cosh (y/x) 



[x + (1 + z)x D ] 2 - Axzx D s ^f J [1 + x(l + z)y D } 2 - 4xzy D 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



with xd — m 2 D Xf and yr> = m 2 D /s. 

As mentioned before, R 3 depends on baryo-chemical potential /x through the Debye 
screening mass which we evaluate using Juttner, BM and MFD distributions. 
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FIG. 2. (a)The gluon production rate, R3/T as function of X g (X q = X q = X g /5) with Juttner, 
Boltzmann and MFD distribution function, for x = 0. (b) The gluon production rate R3/T as a 
function of x for different values of X g where (X q = X q = X g /5). 

Fig. 2(a) shows R3/T as a function of A 9 for a typical gluon rich plasma. At small values 
of X g , R3/T is not much sensitive on various distribution functions, but they start deviating 
when X g exceeds ~ 0.5. Since the plasma remains highly unsaturated by the time T drops 
to T c , results based on MFD or BM distribution will not differ much from the actual result. 
Fig. 3(b) shows the variation of R3/T with x at few typical values of X g . Again, the x 
dependence on R3 is found to be insignificant particularly at small baryon density. 



C. Chemical equilibration 

The space time evolution of the QGP is studied using ideal hydrodynamics in (1+1) 
dimension along with the following master rate equations |3| 

dn„ 



Or 



(i?2^3 — -^3^2) — (Rg^q — Rg->g) (36) 
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-g^ = -gf = ( R 9-*1 ~ R 1^ 9 ) ( 37 ) 

where R2-,?, and ^3^2 denote the rates for the process gg — > ggg and its reverse and R g -> q 
and Rq^g for the process gg — ► gg and its reverse respectively. Similarly, the equation for 
the conservation of energy and momentum in (1+1) dimension can be written as 

|^ + ^ = (38) 
Or r 

which does not include viscosity effect [|13|]. In case of an ideal fluid, 3p =e. From the 
conservation of baryon number, one gets d^nsu^) = which results in 

n B r = {n q — riq)T = const (39) 

We evaluate the RHS of Eq.(36) and Eq.(37) numerically using Juttner distribution. 
However, in order to solve the above coupled equations, we use factorized density distri- 
butions n g = Xgfbg and n q ^ = ^Q(Q)n q (q) (where Hi are the equilibrium parton densities) 
obtained on the basis of Eqs. (3) and (4). Since the plasma is unsaturated, we could have 
also used the BM approximation (Eq. 5). But, we prefer to use Eqs. (3) and (4) for 
consistency with earlier works 

Finally, we solve the following set of equations for X g , Xq, Xq and T numerically by using 
fourth order Runga-Kutta method @], 

^ + f + — ^3(1-^-2^(1-^), (40) 

Xg T T X 2 g 



Xq , 3f 1 „ di X a ^ XqX 

9 



Xq T t 61 X Q Xi 



Ag = A Q + ^ 7 ^ + ^(A Q -A ), (42) 



3T 1 3 • 

jr + - = ~ ^Ms + Mq + Mq], (43) 



where 



A t = a 2 X g + b 2 {X Q + Xq). 

As in Ref. we take the initial conditions as T = 0.57 GeV, A g o=0.09, A g o=0.02 at 
t =0.3 fm/c and treat Xq{ix/T) as a parameter. 
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FIG. 3. (a) The temperature, the quark chemical potential x = fi/T, the gluon fugacity X g , 
the quark fugacity X q , as a function of r for R2 integrated using Juttner distribution and also 
for B.2 using approximate formula of Biro et al. [4] with the initial conditions To = 0.57 GeV, 
X g o = 0.09 and A^o = 0.02, xo = 1.0. (b)The time dependence of T, x, X g and X q with initial 
conditions T = 0.57 GeV, X g0 = 0.09 and X q0 = 0.02 for x = 0.0, 1.0, 2.0. 

Fig. 3(a) shows temperature T and the fugacities as function of r obtained with the 
revised values of R2 and -R3 (solid curves). The dashed curves are obtained using Biro 
et al's factorization. Since R2 values are found less as compared to Biro et al's result, 
the equilibration rate for \ q slows down and the rate for \ g goes up marginally at large 
t. The equilibration rate for x also slows down. This implies that the rate of decrease 
(increase) of quark (anti-quark) contents becomes slower. However, as increase in the rate 
for X g and decrease in the rate of X q have opposite effects on temperature, the variation 
of T with r remains practically unchanged. In Fig. 3(b), we show the effect of baryon 
density on chemical equilibration with revised rates R2 and R3. As before [0], the gluon 
equilibration rate slows down much more than the slight enhancement found in quark (anti- 
quark) equilibration rate. The overall effect is that the temperature of the plasma falls at a 
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slower rate as compared to the plasma when it is baryon free. 



III. THERMAL HARD PHOTON PRODUCTION 

Hard photons are the promising probes to study the evolution of the quark gluon plasma 
produced in relativistic heavy ion collisions. So far photon emission has been considered 
mostly from either chemically equilibrated or non-equilibrated plasma at zero baryon density 
|I^16[. In Ref. |P7yiB|l , hard photon yield has been calculated from a chemically equilibrated 
baryon rich plasma. Here we investigate the photon production rate and integrated yields 
from a chemically unsaturated plasma at finite baryon density. 

In a thermodynamic system, the photon production rate and its momentum distribution 
depend on the momentum distribution of the quarks, anti-quarks and the gluons, which 
is governed by the thermodynamic state of the plasma. To get the emission rate of a 
photon from the plasma we have to fold the amplitude for these reactions with the thermal 
distribution functions and integrate over the phase space volume of all particles except 
photon . The calculation of photon emission rate from a QGP follows via two steps. The 
emission rate from a stationary plasma at a fixed temperature T is determined using Juttner 
distribution function and the yield is calculated by integrating over the plasma volume 
created by the expansion and consequent cooling of the plasma. In thermal and chemical 
equilibrium, the production rate of hard photons with energy E ^> T can be computed 



using the Braaten-Yaan prescription |l9fl. Then the rate decomposes into a soft part, which 



is treated using a resummed quark propagator according to Braaten-Pisarsky Method | 
and a hard part containing only bare propagators and vertices. The medium effects in the 
QGP are included in the soft part through the thermal quark mass, which serves as an 
infrared cutoff in the case of a vanishing bare quark mass. The hard part follows from 
the momentum integration over the matrix elements that lead to photon emission in lowest 



order multiplied by the distribution functions of the incoming and outgoing partons [TH 



A separation scale k c is introduced, which allows one to distinguish between soft and hard 
momenta of the intermediate quark. 

In the following, we only consider the hard part of the photon emission rate and the 
cut-off parameter k 2 is replaced by the thermal mass 2m 2 . The main processes of photon 
production are ( 0(aa s ) annihilation qq — > gj and Compton processes q(q)g — > q(q)j)- The 
thermal rates for these reactions are 

dRhard N d 3 pi ^ ^ 
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*(p? + p£- p£- ^)EI-^| 2 (44) 

where /i,2,3 are the parton distribution function corresponding to two processes [|7|], with 
plus sign for annihilation and the minus sign for the two Compton processes. Following Ref. 
UnyiSfl , the above equation can be written as 

dR hard 1 r°° , r k c 



2E^- — = - \ 7 / ds I " dtY\M\ 2 
d 3 p 8 (2vr) 7 £ yjjfcg J-s+ki ^ 1 1 



x 



dEi (45) 

« 2 ' JP(E U E 2 ) 



where P(E 1 , E 2 ) = aE\ + 6i? 2 + c, with a,b,c as given by 

a = -(a + t) 2 ; b = 2{s + t){Es - E 2 t) ; c = st(s + t) - (E s + E 2 f (46) 

and is the step function, s, t, u are the Mandelstam variables. The matrix elements for 
the Compton and annihilation processes can be written as 



V \M\ 2 = - — 7r 2 aa s U +t 
1 9 ut 

2 9 R s 2 4- 

El^| 2 = -^vr 2 aa s ^l (47 ) 

For the case \i = 0, it is possible to obtain an approximate analytic expression for Eq.(45) 
assuming BM distribution for fi and f 2 and full Fermi/Bose distribution functions for / 3 ||14|| . 
Since we consider baryon rich unsaturated plasma and use Juttner distribution function, we 
have to employ numerical methods. Following [18], we rewrite Eq.(45) in a form suitable for 
Gauss quadrature with E + := E\ + E 2 

J ID hard Krun/ roo 1 r—k 2 



r dE + e -(z + -z-TE)/T . 
Je +te 



Lagauerre Legendre 

l 



+ ' 1 =F A 3 e-( E + - E )/ T 

v ' 

Lagauerre 

dE 2 Aj A 2 , . 

\(E + ,E 2 ) l±Xie~( E + ~ E ^ T l + \ 2 e- E */ T 1 ) 

-v ' 

Chebyshev 

where the upper signs and Ai = Xq, X 2 = Xq , A 3 = A 9 , is used for annihilation process, the 
Compton processes require lower signs and Ai = A 9 , X 2 = A3 = Xq^q-j and the two results 
for Xq and Xq are to be added in order to take quarks and anti-quarks into consideration. 
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The polynomial P 1 (E + ,E 2 ) in Eq.(48) is P(E + — E 2 ,E 2 )/s 2 and has proper weight for 
Gauss- Chebyshev quadrature in E 2 , where E 2 and E 2 are the two roots of the polynomial 
P(E + , E 2 ). The E + integral as well as the s integral, is done numerically by Gauss-Laguerre 
quadrature and the t integral is performed using Gauss-Legendre quadrature. This way, 
the accuracy of the integral has been improved. We have also verified the result using BM 
distribution for f\ and f 2 that matches with the analytic results corresponding to equilibrated 
plasma at x = in the limit k 2 — > 0. 
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FIG. 4. (a) The photon production rate 2EdR/d 3 p for a chemically equilibrated plasma as a 
function of photon energy E at a fixed initial temperature To = 0.265 GeV corresponds to energy 
density e = 9.0 GeV/fm 3 when x = for x = 0.0, 1.0, 1.5 and 2.0 and X g = X q = 1.0. (b) The 
photon production rate 2EdR/ 'd 3 p as a function of photon energy E at a fixed initial energy density 
e = 9.0 GeV/fm 3 for x = 0.0, 1.0, 1.5 and 2.0 and \ g = X q = 1.0. 
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Figs. 4(a) and 4(b) show the photon production rate for a chemically equilibrated plasma 
both at a fixed temperature (T=0.265 GeV) and fixed energy density (e=9.0 GeV/fm 3 ). The 
energy density e of 9.0 Gev/fm 3 corresponds to an initial temperature of 0.265 GeV at x=0. 
These results are similar to what was found in Ref . fl8| . We have further extended the work 
to calculate the photon production from unsaturated plasma at finite baryon density. Figs. 
5(a) and 5(b) shows the corresponding results for an unsaturated plasma. The fugacities 
shown in the figure caption correspond to an initial temperature of 0.57 GeV for energy 
density e=9.0 GeV/fm 3 at x — 0. Both Figs. 4(a) and 5(a) show similar behaviour namely, 
that the rate of photon production is enhanced with increasing baryon density. As will 
be shown below with analytic expression, the annihilation part is not affected much by x 
whereas the Compton process for quark(anti-quark) is enhanced (suppressed) exponentially 
(e ±x ). Therefore, the enhancement with increasing baryon density comes due to the quark 
part of the Compton process. At fixed energy density, the temperature also decreases with 
the increasing baryon density. Hence the net effect is that the photon rate decreases with 
increasing x if the plasma is saturated (see Fig. 4(b). In the case of unsaturated plasma, 
the decrease of initial temperature with increasing x is not very significant to compensate 
the Compton enhancement. The overall effect is, therefore, that the rate still increases with 
increasing x. 

Since we are mainly interested to calculate the integrated photon yield which involves a 
six- dimensional integration (see Eq. 52 below) and will be time consuming, we look for an 
approximate analytic solution for photon production rate i.e. an approximation to Eq.(45). 
For a chemically unsaturated plasma, we can use MFD distribution functions to factorize 
the product appearing in Eq. (44) as 

h / 2 (1 ± h) = AiA 2 A 3 A e<? f? (1 ± f e , q ) + AiA 2 (l - A 3 )/f ft q . (49) 



The factorization is similar to what has been used in Ref. JTBJ to estimate the photon yield 
from a baryon free unsaturated plasma except that the fugacities Ai , A 2 , A3 are replaced 
by X g , \q , Xq corresponding to gluon, quark and anti-quark respectively. The photon 
production rate can now be written in two parts. 
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FIG. 5. (a) The photon production rate 2EdR/d 3 p chemically non-equilibrated plasma as a 
function of photon energy E at a fixed initial temperature To = 0.265 GeV for x = 0.0, 1.0, 1.5 and 
2.0 and \ g = .09, \ q = .02. The dots are the results of the approximate formula (sum of Eq.(51) 
Eq. (52)). (b) The photon production rate 2EdR/d s p as a function of photon energy E at a fixed 
initial energy density e = 9.0 GeV/fm 3 for x = 0.0, 1.0, 1.5 and 2.0 and X g = .09, X q = 1.0. The 
dots are the results of the approximate formula (sum of Eq.(51) and Eq. (52)). 

For the first part, one can use the analytic form that can be obtained using BM dis- 
tribution for /i and fi as in the case for equilibrated baryon free plasma [|T^] and given 
by 



2E 



5a a* 



d 3 p ) x 9 7T 2 



T 2 e~ E l T 



4 ET 



1.43 



annihilation 

'4 ET\ 



k 2 c J 



+ 0.0075 



Compton 



(50) 



Following Ref. flTofl, the second part can be written as 
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' 2E— \ = 10a ° s T 2 e~ E/T 

. d 3 P L ~ 9 7T 4 



X q \q(1-X 9 ) -2-2 7 + 2/n 



'4 ET\ 
7W) 



annihilation 



+ \ {AqA 9 (1 - X Q ) + X Q X g (l - X Q )} |l - 2 7 + 2 In (^^j } (51) 
v , ' 

Compton 

where a, a s are electromagnetic and strong coupling constants respectively. In the above, 
the first term is the contribution from the annihilation whereas the second and third terms 
are due to Compton like processes. The soft part of the contribution is taken care by 
choosing the cut-off parameter k 2 = 2m 2 . As argued before, the annihilation part is nearly 
independent of x (the dependence on /i through k 2 is logarithmic) where as the Compton 
like processes depend on x directly through Xq or Xq. The photon rate is the sum of Eq.(50) 
and Eq.(51). This formula reproduces the results of photon production rate quite accurately 
(when the plasma is unsaturated) with the exact result obtained from numerical integration. 
The results are shown in Figs. 5(a) and 5(b) (see solid circles) for comparison at two values 
of x = and x = 2.0. The above analytic expression also reproduces the result for saturated 
plasma although the agreement is not so exact as in case of Figs. 5(a) and 5(b). 

Next we calculate the photon yield by integrating the above analytic expression of photon 
production rate (sum of Eq. (50) and Eq. (51)) over the space and time (TB| 



/ 2 dn \ r T c rv^c / dn \ 

\dPpTdy~ ) =Q L Lnuc V [ 2E d^ ) El ° C reSt=P - ^y-y'^=T(r),X=X(r) (52) 

The r dependence of T(r), A 9 (r) and Xq^(t) are taken from the above calculations (in 
previous section) with a 5 =0.3, X g0 =0.09, X q0 =0.02 and T = 0.57 GeV. Further, we take 
y nuc = ±6 corresponding to RHIC energy. The transverse cross-section of the gold nuclei, 
Q w 180fm 2 . Figs. 6(a) and 6(b) show the above integrated yields both at a fixed initial 
temperature and fixed energy density. The integrated yields both at fixed temperature and 
fixed energy density behave similarly as the rates shown in Figs. 5(a) and 5(b). It may be 
pointed out that, the finite baryon density effects which result in a slower cooling of the 
plasma affect the thermal photon and dilepton production yields somewhat differently. For 
example the basic rate of thermal dilepton production (studied in previous work Ref. [0) 
does not depend on x as it involves only XqXq. 
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FIG. 6. (a) The integrated photon yield dnjd p T dy from a chemically non-equilibrated plasma 
at a fixed initial temperature To = 0.57 GeV for various values of xq. (b) The integrated photon 
yield from a chemically non-equilibrated plasma at a fixed initial energy density eo = 9.0 GeV/fm 3 
for various values of xq. 



But the integrated yield is enhanced due to the slower space time evolution at finite 
baryon density. On the other hand, in case of thermal photon production, the Compton like 
processes involve the product of the type \ 9 \q or A 9 Aq. Therefore, unlike the dilepton case, 
the basic photon yield also gets affected due to chemical equilibration since the variation of 
x and \ g are significant. 



IV. SUMMARY AND CONCLUSION 



In the present work, we have investigated the effect of baryon density on parton pro- 
duction for the process gg # ggg and gg # qq with full phase space distribution and also 
including quantum effects like Pauli blocking and Bose enhancement factors explicitly. The 
results obtained using exact distribution function (Juttner) does not differ much from the 



21 



results obtained using factorized distribution function of MFD or Boltzmann (BM) type. 
Although, BM approximation is classical, the final expression for the rate gg # qq includes 
quantum statistics. The results which are obtained on the basis of Biro et al's estimation, 
although classical in nature, overestimate the rate by 20% to 30% due to use of MFD type 
of approximation. Further it is found that, the rates for both the processes {gg # ggg and 
gg # qq) are insensitive to baryon density particularly when the plasma is unsaturated. 
This further justifies the use of MFD or BM type of distribution for partons for unsaturated 
plasma at small baryon density. We have also studied the process of chemical equilibration 
using the revised rates and also compare the results with the previous values obtained as 
per Biro et al's estimation. It is found, that the quark and gluon equilibration rates are 
affected slightly whereas the temperature remains unaffected. We have also studied thermal 
hard photon yields from a chemically non-equilibrated plasma which are considered to be 
the ideal probe to study the space time evolution of the plasma. We employ numerical tech- 
niques to estimate photon rate using Juttner distribution functions and also compare our 
results with approximate analytic expression. It is seen that as in the case of the dilepton 
yields, the photon yield is also enhanced in presence of finite baryon density. 
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